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We utilise a form for the Hubble parameter to generate a number of solutions to the Einstein field 
equations with variable cosmological constant and variable gravitational constant in the presence of a 
bulk viscous fluid. The Hubble law utilised yields a constant value for the deceleration parameter. A 
new class of solutions is presented in the Robertson- Walker spacetimes. The coefficient of bulk vis- 
cosity is assumed to be a power function of the mass density. For a class of solutions, the deceleration 
parameter is negative which is consistent with the supernovae la observations. 

1. Introduction 

Einstein proposed his General Theory of Relativity as a geometric theory. He introduced 
a cosmological constant into his field equations in order to obtain a static cosmological 
model as without the cosmological term, the field equations would admit only non-static 
solutions. The Einstein field equations are a coupled system of highly nonlinear differential 
equations and we seek physical solutions for applications in cosmology and astrophysics. 
In order to solve the field equations we normally assume aform for the matter content or 
suppose that spacetime admits Killing vector symmetries El. Solutions to the field equa- 
tions may also be generated by applying a law of variation for Hubble's parameter. It is 
interesting to observe that this law yields a constant value for the deceleration parameter. 
The variation of Hubble's law assumed is not inconsistent with observation and has the 
advantage of providing simple functional forms of the scale factor. In the simplest case 
the Hubble law yields a constant value for the deceleration parameter. In earlier literature 
cosmological models with a constant deceleration parameter have been studied by Berman 
!■!, Berman and Gomide i, Johri and Desikan i, Singh and Desikan i, Pradhan et al and 
others. 



Models with a relic cosmological constant A have received considerable attention re- 
cently among researchers for various reasons (see Refs. and references therein). Some 
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of the recent discussions on the cosmological constant "problem" and on cosmology with a 
time-varyingcosmological constant by Ratra and Peebles 113, Dolgov Ilia Il3 and Sahni and 
Starobinsky E_ZJ point out that in the absence of any interaction with matter or radiation, the 
cosmological constant remains a "constant", however, in the presence of interactions with 
matter or radiation, a solution of Einstein equations and the assumed equation of covari- 
ant conservation of stress-energy with a time-varying A can be found. For these solutions, 
conservation of energy requires decrease in the energy density of the vacuum component 
to be compensated by a corresponding increase in the energy density of matter or radiation. 
Earlier researchers on this topic, are contained in ZeldovichEJ, Weinberg El and Carroll, 
Press and Turner 0. Recent observations by Perlmutter et al t§ and Riess et al strongly 
favour a significant and positive A. Their finding arise from the study of more than 50 type 
la supernovae with redshifts in the range 0.10 < z < 0.83 and suggest Friedmann models 
with negative pressure matter such as_a cosmological constant, domain walls or cosmic 
strings (Vilenkin ca, Garnavich et al cj. The main conclusion of these works is that the 
expansion of the universe is accelerating. 

SeveraLansatz have beenjiroposed in which Jhe A term decays witLlime (see Refs. 
Gasperini E-aEa, Freese et al c3, Qz.er and Taha 113, Peebles and Ratra 123, Chen_and Hu 
c3, Abdussattar and Viswakarma c3, Gariel and Le Denmat E2l, Pradhan et al til ). Of 
the special interest is the ansatz^i oc S~ 2 (where S is the scale factor of the Robertson- 
Walker metric) by Chen and Wu c3, which has been considered/modified by several authors 
( Abdel-RahamanE3, Carvalho et al L3, WagaH, Silveira and WagaE-3, VishwakarmaE3). 

In most treatments of cosmology, cosmic fluid is considered as perfect fluid. However, 

Slk_yiscosity is expected to play an important role at certain stages of expanding universe 
It has been shown that bulk viscosity leads to inflationary like solution Lj, and acts 
like a negative energy field in an expanding universe 113. A number of authors have dis- 
cussed cosmological solutions with bulk viscosity in various context 

It has been shown by Berman i and Berman & Gomide i that all the phases of the 
universe, i.e., radiation, inflation and pressure-free, may be considered as particular cases 
of the deceleration parameter q = constant type, as 

SS 

q = -^ 

where dots stand for time derivatives. We extend this definition to the Robertson- Walker 
cosmological models. In the past few decades there have been numerous modification of 
general relativity in which gravitational constant G varies with time llj. Considering the 
principle of absolute quark confinement, Der Sarkissian @ has suggested that gravitational 
and cosmological constants may be considered as functions of time parameter in Einstein's 
theory of relativity. A number of authors have considered time-varying G and A 

within the frame work of general relativity. 
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Maharaj and Naidoo E3 utilised a form of the Hubble parameter to generate a number of 
solutions to the Einstein field equations with variable cosmological constant and variable 
gravitational constant in presence of a perfect fluid as the source of matter. These authors 
have extended the results obtained by Berman H, Berman and Gomide □ by obtaining so- 
lutions to Einstein field equations, with variable gravitational and cosmological constants, 
in the Robertson-Walker spacetime. Explicit forms for the gravitational constant, cosmo- 
logical constant, scale factor, energy density and pressure are obtained for various cases. 



Motivated by the situations discussed above that bulk viscosity, gravitational and cos- 
mological "constants", are more relevant during early stages of the universe, our intension 
in this paper is to extend the results obtained by Berman H, Maharaj and Naidoo EJ by 
including a bulk viscous fluid as a source of matter in the energy momentum tensor. This 
paper is organized as follows. In section 2, the solutions and main results of Berman i, Ma- 
haraj and Naidoo E3 are reviewed. In section 3, the corresponding solutions for a universe 
filled with a bulk viscous fluid are found. In section 4, a number of classes of new solutions 
for all cases of k : 0,1,-1 for variable G and A are presented. Our results are discussed in 
section 5. 



2. Robertson- Walker Spacetime-Revisited 

In the standard coordinates (x a ) = (t, r, 8, 4>) the Robertson- Walker line element has 
the form 

dr 2 



ds 2 



-dt 2 + S 2 (t) 



1 — kr 2 



+ r 2 (d9 2 + sin 2 9d<j) 2 ) 



(1) 



where S(t) is a cosmic scale factor. Without loss of generality the constant k is related to the 
spatial geometry of a 3-dimensional manifold generated by t = constant. The Robertson- 
Walker spacetimes are the standard cosmological models and are consistent with obser- 
vational results. For the case of variable cosmological constant A(t) and gravitational 
constant G(t) the Einstein field equations 



Gab + Ag ab = 8nGT ab 



(2) 



yield 



3 



{S 2 + k) = 8wGfi + A 



S 2 

.S (S 2 + k) 



s 



s 2 



= SirGp + A 



(3) 



(4) 



for the line element (1). From Eqs. (|3j) and (|4J) we obtain the generalised continuity equa- 
tion 



. . G A 

M + 3-( Ai + p) + - M + — 



0. 



(5) 
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This reduces to the conventional continuity equation for constant A and G. If the clas- 
sical conservation law, T a ^ = 0, also holds, then Eq. (||) implies two relationships 

/i + 3-(/x+p) = 0, (6) 

8-khG + A = 0, (7) 

which facilitate the solutions of the field equations. The result (^) is just the conventional 
continuity equation, and (Q) simply relates G and A and does not explicitly contain the 
scale factor S(t). 

Maharaj and Naidoo considered the generalised Einstein field equations (||)-(0)with 
variable gravitational constant G and variable cosmological constant A (i) for the Robertson- 
Walker metric (1). They assumed the variation of the Hubble parameter as given by the 
equation 

H = - = DS- m (8) 

where D and m are constants. Then this imply that the deceleration parameter q is constant 
i. e. q = m — 1. 

The form of the Hubble parameter was first utilised by Berman and Berman and 
Gomidea for the case of classical Einstein field equations with A = and G = 0. Berman 
u presented a solution to the field equations (Q) for the k = Robertson- Walker spacetime: 

S= (C + mDt) 1 /™, (9) 

A = BS- 2m , (10) 
G = ps mB K** A \ (11) 

fl= ^ S -2m-mB/^A)^ (u) 
P 



A 



B 

m(2 



g — 2m — 7715/(471-^4) (13) 



An A 

where A, B, C, p are constants and are subject to the following condition 

3D 2 = 8ttA + B 



It isworth noting that Eq. ( |l6| ) given by Berman u corresponding to Maharaj and 
Naidoo Ell equation (pj|), has an incorrect coefficient on the right hand side. Equations (^)- 
dl3| ) comprise the general solution to the generalised Einstein field equations @)-(Q) with 
variable G and A for the Hubble law (||). 
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3. Bulk Viscous Solutions of the Field Equations 



In this section bulk viscous models of the universe are discussed. Weinberge a has 
suggested that in order to consider the effect of bulk viscosity, the perfect fluid pressure 
should be replaced by the effective pressure p by 



(14) 



where £ is the coefficient of bulk viscosity and 9 is the expansion scalar given by 9 = 3H. 
Here £ is, in general, a function of time. 
Therefore, from Eq. (O), we obtain 



m(2 



B 
AttA 



)-3 



C<— 2m— mB/(4,nA) 



(15) 



For complete determinacy of the system, we assume an equation of state of an ideal gas 
given by 

P = 7M, 0< 7 <1, (16) 

where 7 is a constant. Thus, given £(f) we can solve for the cosmological parameters. It is 
standard to assume *M1 the following widely accepted ad hoc law 



(17) 



If n = 1, Eq. ( |17| ) may correspond to a radiating fluid, whereas n = 3/2 may correspond 
to a string-dominated universe H. However, more realistic models E3 are based on ri lying 



(18) 



the regime < n < 1/2. On using Eq. (Ill) in Eq. (1151), we obtain 



A 



m(2 



4ttA 



ct—2m—mB/(4:irA) 



3.1. ModeZ 7; (£ = ft) 

When n = 0, Eq. ( |l7| ) reduces to £ = ft = constant and hence Eq. ( |l8| ) with Eq. ([l6]) gives 



3^ft 



.4 

3/3 



m(2- 



47rA' 



Q—2m—mB/(AirA) 



(19) 



3.2. ModeZ 71; (£ = ft/x) 

When n = 1, Eq. ( |l7| ) reduces to £ = ft/i and hence Eq. ( |l8| ) with Eq. ( Jl6] ) gives 



.4 



3/3(7 - 3Z?ft5-» l ) 



m(2 



B 
A-kA 



)-3 



g— 2rn— mB/(iirA) 



(20) 



It is possible to avoid the horizon and monopole problem with the above variable G(t) and 
A(£l solutions as suggested by Berman B. Other models are also considered by Berman et 
al E3 and Bertolami cZl which have the relationship 



Aoc-. 
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This form of A is physically reasonable as observations suggest that A is very small 
in the present universe. A decreasing functional form permits A to be large in the early 
universe. A partial list of cosmological models inwhich the gravitational constant G is 
decreasing function of time are contained in Gr0n cM, Hellings et al E3, Rowan-Robinson 
ELD, Shapiro et al E3 and Van FlandernEll. The possibility of G increasing with time, at least 
in some stages of thedevelopmentof the universe, has been investigated by Abdel-Rahman 
l3, Chow Ej, Levitt EH and Milne Ex In these models, by selecting proper signs of different 
constants, one can make fi positive and decreasing function of time t. 



4. Other Solutions-Revisited 

Maharaj and Naidoo presented a number of classes of new solutions for all classes of 
k : 0, 1, —1 for variable cosmological constant A and variable gravitational constant G for 
Hubble law (^|). These solutions covered both the cases of to — and ra/0 for the scale 
factor S: 

[C + mDt] 1 /™ whenm^O 
Ee Dt when m = 

where C, D, E are constants. To solve the Einstein field equations (§)-(§) they adopted the 
ansatz 

3D 2 

' K, (22) 



S(t) = 



(21) 



g2m 

8irGfi 



A 
3k 

s 2 



= K 



(23) 



where K is constant. Thisjansatz has the advantage of providing further classes of solu- 
tions. Maharaj and Naidoo Ell obtained the expressions for cosmological constant A, energy 
density /i in terms of the gravitational constant G and the scale factor S 



A = K 



S 2r. 



8irG 



3k „ 



(24) 



(25) 



On substituting Eq. (|25|) and the derivative of Eq. (Q) with respect to time coordinate 
t into Eq. (^|) we obtain the differential equation 



— = 6mD — - 

G S 2 "^ 



1 / 3k 



K) 



(26) 



relating G to S. Thus the scalar factor S is specified by our assumed form of the Hubble 



parameter, the gravitational constant G is known in principle. The ansatz (|22|)-([23[) enables 
to integrate all the Einstein field equations for a number of values of to, k and K. In the 
remainder of this section we present a variety of classes of solutions to the Einstein field 
equations for each of these cases considered in the presence of bulk viscosity. We list the 
form of the scale factor S, the variable cosmological constant A, the variable gravitational 
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constant G, the energy density /i and the effective pressure p. There are other classes of 
solution possible for other values of m. However the integration becomes extremely com- 
plicated for general m, we present only some simple cases in the following. 



4.1. Case J:m = 0,iC^0,^0. 
In this case we obtain 



S = Ee D \ 
A = 3D 2 - K, 
G = A 1 
1 



/ <- : 



8ttA 
1 



3k „ 



k 

S 2 



K 



(27) 
(28) 
(29) 

(30) 
(31) 



4.1.1. Model I :(£ = £o) 

When n = 0, Eq. ( JTv| ) reduces to £ = £o = constant and hence Eq. ( |3~l| ) with Eq. (|l6|) gives 

1 



1 

M = - 

7 



8ttA 



(32) 



In this model, if we set A < and £ ,D,K,k > then /i is always positive and de- 
creasing function with time. 



4.1.2. Model II : (£ = £ /i) 

When n = 1, Eq. ( |l7| ) reduces to ^ = and hence Eq.(|3l]) with Eq. (|l6| gives 

1 

M_ ~87rA( 7 -3i}£o) 

In these de Sitter-type solutions A and G are strictly constant because of the restric- 
tion m = 0. The cosmological constant A vanishes when K = 3D 2 and is positive for 
K < 3D 2 . The scale factor S is exponential in t, so that if D > then the universe is 
exponentially expanding always. Such models are not physical description of our present 
universe but could be applicable in the early universe in the inflationary scenario. For 
m = 0, we get the deceleration parameter q = — 1 for these class of solutions which is con- 
sistent with the recent observations of supernovae la which require that the present universe 
is accelerating 00. In model II, if we set A < 0, and S, ,D,K,k > and 7 > 3D£ , 
then fi is always positive and decreasing function with time. 



A: 
52 



K 



(33) 
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4.2. Case II : m ^ 0, K = 0, k ^ 0. 
In this case we obtain 



S = [C + mDt] 1/T 
3D 2 



A 



S 2r 



k(l — m) 
8na k(l ~ m) 



1 

8ira 



AmD 3 3k 



exp{-^-S 2 ~ 2 ™}. 
k(l — m) 



(34) 
(35) 

(36) 
(37) 
(38) 



4.2.1. Model I : (£ = £ ) 

When 77 = 0, Eq. (|l7|) reduces to £ = £0 = constant and hence Eq. ( |38| ) with Eq. (|l6|) gives 



1 

M = - 

7 



8"7ra 



S 3 ™ 



3fc 
52 



fc(l — 772 J 



(39) 



4.2.2. Mode? // : (£ = £ o a0 

When 77 = 1, Eq. ( |l7| ) reduces to £ = £oM an d hence Eq.(|3~8|) with Eq. ( JT6| ) gives 



1 



87ra( 7 - SZ^oS"-" 1 ) 



4mL> 3 



3k 
^2 



777 D 



fc(i — m) 



S 2 



(40) 



This case shares the common feature that G may be increasing in time in certain regions of 
spacetime with the model proposed by Abdel-Rahman EH. 



4.3. Case III : m ^ 0, K ^ 0, k = 0. 
In this case we obtain 



S =[C + mDt] 
3D 2 



l/m 



A = 



S 2r ' 



P = 



I' = 
1 

8ira 



G = a exp{ 
exp{ 



-K, 
3D 2 

KS 2m 

3D 2 



h 



8na 
2mD 2 
S 2m 



KS 2 



3D 2 
KS 



2)ii 



}■ 



(41) 
(42) 

(43) 

(44) 

(45) 
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4.3.1. Model I : (£ = f„) 

When n = 0, Eq. (|l7[) reduces to £ = £o = constant and hence Eq. (p3h with Eq. (n~6h gives 



1 

M = - 

7 



3D6,S -m - J- 
07ra 



S 2 ™ 



3D 2 



(46) 



4.3.2. Mode/ JJ : (£ = £ m) 



When n = 1, Eq. (|17[) reduces to £ = £oM an d hence Eq. (|45|) with Eq. ( 16) gives 



1 



/' 



87ra(7 - iD^S- 



2mD 2 
S 2m 



K 



-3D 2 



(47) 



4.4. Case IV : m = 2, K ^ 0, k f 0. 
In this case we obtain 

S =[C + 2Dt] 1 ' 2 , 
A = -^r- - if, 



,S 4 

(3fc5- 2 + K ) R / k 



2D 2 /k 



H = \3kS~ 2 + K] 

8na 



p=- — {3kS- 2 + K] 
Una 



exp{S 2 } 

exp{S- 2 } 



(3kS- 2 + K) K l k 
exp{S- 2 } 



4D 2 S- i (2K - 3kS~ 2 ) 



1 



(3kS~ 2 

kS' 2 
Ana 



2D 2 /k 



2D 2 /k 



exp{S~ 2 } 



-, 2D 2 /k 



(3kS~ 2 + K) K l k 



(48) 
(49) 

(50) 
(51) 



(52) 



Unlike the cases considered thus far we have a specific value for m. This gives a value 
q = 1 for the deceleration parameter. A wide range of behaviour is possible for the gravi- 
tational constant. 



4.4.1. Model I : (£ = f ) 

When n = 0, Eq. ( JTv| ) reduces to £ = £o = constant and hence Eq. (|52|) with Eq. (|l6|) gives 

2D 2 /k 



-[3kS- 2 + k] 



7 87ra7 

AD 2 S~ 4 {2K- 3kS~ 2 ) 
k(3kS- 2 + K) 



exp{S~ 2 } 



(3kS- 2 + K) K l k 
2kS- 2 



1 - 



(3kS- 2 + K) 



(53) 
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4.4.2. Model II : (£ = £ o a0 

When n = 1, Eq. ( |l7| ) reduces to £ = £o/i and hence Eq. ( |52| ) with Eq. ( |l6| ) gives 

2D 2 /fc 



1 



_ k{3kS~ 2 + K) 



exp{S~ 2 } 



(3kS- 2 + K) K / k 

2kS~ 2 
~ (3kS~ 2 + K) 



(54) 



4.5. Case V : m = -2, K f 0, k ^ 0. 
In this case we obtain 

S = 



A = -tt—t - K, 



G = a- 



8- 
exp{ 



MS' 
T 2 " 



six 

IK > 



[S-^kS^ + K)} 9 * 2 ^ 3 ' 

a-l&k 2 /K a o4 OJ.C2 



87ra 

p = [3kS- 2 + K] 9k ■ l K + 1 



4na 



6k 2 



2K K 2 
3k -S 6 



K 3 3S 2 (3kS~ 2 + K) 2 



6XP{ 2K ~ 



(55) 
(56) 
(57) 

(58) 



(59) 



4.5.1. Model J: (£ = £ ) 

When n = 0, Eq. ( |l7| ) reduces to £ = £o = constant and hence Eq. (^) with Eq. (|l6|) gives 

3^ ^ 2 ^- 18fc2 /^ 3 



/' 



6k 2 



3k -S 6 



[3kS~ 



1 



K 3 3S 2 (3kS- 2 + K) 2 



K 



5' 4 



2^ K 2 



}• 



(60) 



4.5.2. Mode? II : (£ = 

When n = 1, Eq ( |l7| ) reduces to £ = £ M an d hence Eq. ( |5^ ) with Eq. ( |l6| ) gives 

g-i8fc 2 /.ff 3 



M = 



6fc 2 



47ra(7 - 3DS, S 2 ) 



3k -S 6 



K 3 3S 2 (3kS- 2 + K) 



[3kS~ 

1 

2 



(61) 
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For m = —2, we get the deceleration parameter q = —3 for these class of solutions which 
is consistent with therecent observations of supernovae la which require the present uni- 
verse is accelerating E30 



4.6. Case VI : m=\,K^ 0, k ^ 0. 

1 



2 

In this case we obtain 

' 12 



S= [C+-Dt] 2 , (62) 
A = K, (63) 



G = aexp{^-^D 2 arctan( ] J^-S)} 7 (64) 



»=T~ fi kS ~ 2 + K ] ™P{-\ ^-D 2 arctanU §-S)}, (65) 
8na V kK V ok 

p = - J- [(k - D 2 S)S- 2 + K] exp{-\l^-D 2 arctanuf^-S)}. (66) 
57ra V kK V 3fc 

For this class of solutions the deceleration parameter has the value q = asm = ~. 



4.6.1. Model I : (£ = £ ) 

When n = 0, Eq. ( JTt| ) reduces to £ = £o = constant and hence Eq. ( ]66| ) with Eq. (|l6|) gives 

/ 3 n<[ 

exp{-\l —D 2 arctan{\l —S)}. (67) 

4.6.2. Model II : (£ = £ ^) 

When n = 1, Eq. ( |l7| ) reduces to £ = £o/^ ar, d hence Eq. (|66|)with Eq. (|l6|) gives 

exp{-J^D 2 aTctan(J^S)}. (68) 
For these class of solutions, it is observed that the universe is accelerating. 



4.7. Case VII : m = §, K + 0, k + 
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In this case we obtain 



G = 



S = [C+^Dt] 3 ' 2 , 

exp{6 arctan{ -= )} 



£4/3 _ aS 2IZ + a 



H= - — [3kS~ 2 + K] 



£4/3 _ aS 2/3 + a 



1 

87TQ! 



(S 2 / 3 + a) 2 

£4/3 _ a g2/3 + a 2 



V3o 



2S 2/3 - a ' 
exp{— 6 arctan{ — )} 



V3< 



a 



(69) 
(70) 

(71) 

2 

I 

) 

(72) 



r o ,2s 2 / 3 -a,.' 

(g2/ 3 + Q)2 e ^i-6 ar C tan(— ^— )} 

(S 2 / 3 + a) 2 " 
54/3 _ a £2/3 + a 2 



2D 2 

kS- 2 + K+ ——(3kS- 2 + K) 
6Ka 



S 2 ' 3 



4^3(5 4 / 3 - aS 2 ' 3 



(S 2 / 3 + a) 3 (S 2 / 3 + a) 2 {3a 2 + (2S 2 / 3 - a) 2 } 



(73) 



Here the deceleration parameter has the value q = — | as m = | . For these class of solu- 
tions, it is observed that the universe is accelerating. 



4.7.1. Model I : (£ = £ ) 

When n = 0, Eq. < JTv| ) reduces to £ = £0 = constant and hence Eq. ( f73| ) with Eq. (|l6|) gives 



3ggo 1 
7S 2 / 3 8iraj 



£4/3 _ a5 2/3 + ffl 



(S 2 / 3 + a) 2 



2 r , ^S 2 ' 3 -a^ 
exp{— 6 arctan{ r= )} 



2D 2 

fcS" 2 + K + — — (3fc5- 2 + A') 



£2/3 



\/3a 
(S 2 / 3 + a) 2 

£4/3 _ aS 2/3 + a 2 



4V3(5 4 / 3 -a^ 2 / 3 + a 2 ) 

(52/3 + a )3 (52/3 + a )2{ 3fl 2 + ( 2 £2/3 _ a )2} 



(74) 



4.7.2. Model II : (£ = £ o a0 

When n = 0, Eq. |l7|) reduces to £ = £oM an d hence Eq. ( |73"| ) with Eq. (|l6|) gives 



1 



/i = 



54/3 _ a5 2/3 + , 



87ra( 7 - 3D£ S~ 2 / 3 ) [ (S 2 / 3 + a) 2 
2D 2 



-exp{—6 arctan( 



2S 2 / 3 - a 
\/3a 



)} 



D 2 /Ka 2 



3Ka 



(3kS~ z + K) 



(S 



2/3 



£4/3 _ aS 2/3 + a 2 
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S 2 ' 3 - a _ 4V3(5 4 / 3 - aS 2 ' 3 + a 2 ) 
(S 2 / 3 + a) 2 (S* 2 /3 + a) 2 {3a 2 + (2S 2 / 3 - a) 2 } ' (? } 

In the above we have presented a number of new solutions to the Einstein field equa- 
tions with variable cosmological constant and gravitational constant which satisfy the Hub- 
ble variation law given by Eq. (Q). It remarkable that this simple law leads to a wide class of 
solutions. It is interesting to observe that solutions are admitted in which the gravitational 
constant may be increasing with time (cf. Abdel-Rahman, 1990). The anstaz utilised to 
solve the Einstein field equations (3)-(4) is very simple. It might be worthwhile to investi- 
gate other possibilities that lead to solutions to the Einstein's field equations with interesting 
behaviour of the gravitational constant and cosmological constant. 



5. Conclusions 

In this paper we have investigated Einstein's equation in the presence of a viscous fluid, 
for the Robertson-Walker universe within the framework of general relativity, where the 
gravitational constant G and the cosmological parameter A are variables. We utilize a form 
for the Hubble parameter (H = DS~ m ) to generate a number of solutions to Einstein field 
equations. For these class of solutions where m = 0, —2, |, |, we find the decelerating 
parameters as negative. These class of solutions are consistent with the recent observations 
of supernovae la EJEII which require the present universe is accelerating. From our results 
we observe that fi is decreasing with time with suitable choice of constants whereas G 
is increasing function of time. The possibility of an increasing G has been suggested by 
Abdel-Rahaman E3, Arbab EJ and Massa E± 

Assuming an ad hoc law of the form £(t) = £ofJ. n , where /i is the energy density and 
n is the positive index, we have obtained exact solutions. The models discussed here are 
isotropic and homogeneous and, in view of the assumption of isotropy the sheer viscosity 
is absent. The effect of the bulk viscosity is to produce a change in the perfect fluid. We 
observe that Murphy's conclusion E3 about the absence of big bang type singularity in the 
finite past in models with bulk viscous fluid is, in general, not true. 
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